Detecting a small perturbation through its non-Gaussianity 
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A highly non-gaussian cosmological perturbation with a flat spectrum has unusual stochastic 
properties. We show that they depend on the size of the box within which the perturbation is 
defined, but that for a typical observer the parameters defining the perturbation 'run' to compensate 
for any change in the box size. Focusing on the primordial curvature perturbation, we show that an 
un-correlated gaussian-squared component is bounded at around the 10% level by the WMAP bound 
on the bispectrum, and we show that a competitive bound may follow from the trispectrum when 
it too is bounded by WMAP. Similar considerations apply to a highly non-gaussian isocurvature 
perturbation. 

PACS numbers: 98.80.Cq 



Introduction. The origin of structure in the Uni- 
verse seems to be the primordial curvature perturbation 
(, present already a few Hubble times before cosmologi- 
cal scales enter the horizon and come into causal contact 
0. Within the observational uncertainty, C is Gaussian 
with a practically scale-independent spectrum. Future 
observation though may find non-Gaussianity or scale- 
dependence and in this Letter we focus on the former. 
We shall discount the possible contribution of topological 
defects to observables, but we shall consider briefly the 
possible contribution of a primordial isocurvature pertur- 
bation. 

We assume that £ originates from the vacuum fluctua- 
tions during slow-roll inflation of one or more light scalar 
fields. These fluctuations are promoted to practically 
gaussian classical perturbations around the time of 
horizon exit. Expanding £ in powers of the field perturba- 
tions, the linear terms are Gaussian and quadratic terms 
are expected to account adequately for non-Gaussianity 



To obtain detectable non-Gaussianity some field other 
than the inflaton must contribute to f 0. 0] Possible 
examples of such a field in the literature are the extra 
field in double inflation [i^, the curvaton (see also 
0), and a field causing inhomogeneous reheating Q or 
preheating |2^ . 

When discussing non-Gaussianity, it is usually sup- 
posed that the contribution of the extra field will actually 
dominate £. Then £ is of the form 

C(x)=a(x)-^/ NL (a 2 (x)-( ( x 2 (x))) , (1) 

where cr(x) is the perturbation of the extra field, normal- 
ized so that the linear term has unit coefficient. 1 



We will suppose instead that the additional contribu- 
tion is sub-dominant so that 

C(x) = Cg(x) + C.(x) (2) 
C(x) = aa(x)+a 2 (x)-(a 2 ), (3) 

where Cg is the practically Gaussian contribution of the 
inflaton, and we now normalize a so that the quadratic 
term has unit coefficient. For the moment we take the 
linear term to be negligible so that 2 

C CT (x)=a 2 (x)-(a 2 ). (4) 

Basic Definitions. To describe the stochastic proper- 
ties of cosmological perturbations one formally invokes 
an ensemble of universes, of which the observable Uni- 
verse is supposed to be a typical member. A sampling 
of the ensemble may be regarded as a sampling of dif- 
ferent locations for the region under consideration. The 
stochastic properties are conveniently described using a 
Fourier expansion, which we make inside a finite box of 
coordinate size L much bigger than the region of inter- 
est. In terms of physically significant wavenumbers, this 
means k 3> L^ 1 . 

We will denote a generic cosmological perturbation, 
evaluated at some instant, by g(x), and assume (g) = 0. 
The fields responsible for £ are smoothed on some scale 
^max an d we consider the era when this scale is outside 
the horizon. Taking the box size to be formally infinite, 
the spectrum P g (k) is defined by 

<5kffk'> = (2^ 3 >(k + k')P 3 (fc), (5) 

where gk = / d 3 x e lk x </(x) are the Fourier modes of 
,g(x). The S- function and the dependence of the spec- 
trum only on k = |k| express the fact that stochastic 



1 We are defining the curvature perturbation as inj|| , and the fac- 
tor —3/5 appears because the original authors 9] worked (using 
first-order perturbation theory) with $ = — §C- 



2 The literature includes examples where the linear term dominates 
llll and where the quadratic term dom inates ll J . The most 
closely related previous discussion is ll.'j , whose authors consider 
Eq. Q but with a sharply peaked spectrum for a instead of the 
usual flat spectrum. 
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properties are invariant under rotations and translations. 
The variance is 



<S 2 (x)> 



1 



d 3 kP g {k) 



dk 

T 



Vgik), (6) 
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and V g = (k 3 /2Tr 2 )P g is the typical value of g 2 . On cos- 
mological scales Vq is almost scale-invariant with V^ 2 — 
5 x 10~ 5 . 

Non-Gaussianity is signaled by additional connected 
correlators. The bispectrum B g is defined by 

{Sk^k.Ska} = (2^) 3 ^ 3 )(k 1 +k 2 +k 3 )B s (fc 1 ,A :2 ,fc3)- (7) 

If the curvature perturbation has the form its bispec- 
trum is given to leading order by [9j 

B c {k x ,k 2 M) = -~/nl [P C (ki)Pc(k2) + cyclic] . (8) 

(Only the term linear in /nl is kept, which is justified 
because the second term of Eq. is much less than the 
first term.) Current observation 0,0] gives |/nl| ^ 
100, which makes the non-Gaussian fraction of £ less than 
100P C 1/2 - 10~ 3 . Absent a detection, PLANCK [HI will 
bring this down to roughly |/nl| < 1. 

Following Maldacena Q, we define /nl by Eq. (jHJ ir- 
respectively of its origin, making it in general a function 
of the wavenumbers. 

The trispectrum T g is defined in terms of the connected 
four-point correlator by as 

(ffk^k^k.,^ = (2n) 3 6^(k 1 + k 2 + k 3 + k 4 )T g . (9) 

It is a function of six scalars, defining the quadrilateral 
formed by {ki, k 2 , k 3 , k4}. If the curvature perturbation 
is given by , its trispectrum to leading order is of the 
form |15j 

T c = ir NL [P < (fc 1 )P c (fc 2 )P c (|fc 14 |) + 23perms.] , (10) 

with tnl = 36/^l/25. In this expression, kij = |kj + 
kj|, and the permutations are of {ki, k 2 , k 3 , k,j} giving 
actually 12 distinct terms. 

The authors of 0, H(| have measured some connected 
terms of the angular trispectrum on the COBE DMR 
data, with no significant detection of signals. Non- 
detection of the trispectrum by the COBE DMR data 
may be interpreted as |/nl| < 10 4 (95% c.l.) based upon 
theoretical expectations (Figure 3 of [15j with £ max = 10). 
Expressed as a bound on tnl, the current bound be- 
comes |tnl| ;$ 10 s . When a bound on the trispectrum 
becomes available from WMAP data, it might constrain 
non-gaussianity of the type |JS} more strongly than the 
bispectrum. There is currently no estimate of the bound 
on the trispectrum which will eventually be possible. 

Following Maldacena's strategy for /nl, we will define 
Tnl by Eq. (|10fl irrespectively of its origin, so that in 



general it depends on the wave- vectors. In this Letter we 
will calculate /nl and tnl in the case that the curvature 
perturbation has the form (@J. 

The bispectrum. We take the spectrum V a to be scale- 
invariant, which will be a good approximation if a is 
sufficiently light. Since £ g and a are uncorrelated and the 
curvature perturbation is dominated by C g , the spectrum 
and the bispectrum of the curvature perturbation are 

P c (k) = P Cs (k) + P a2 (k) c P (g (k) (11) 
B ( (k 1 ,k 2 ,k 3 ) = B a 2(k 1 ,k 2 ,k 3 ) (12) 

The Fourier components of a 2 are given by 
1 



(,j2)k = J2nf 
For non vanishing k and k' 



ger q CTk-q • 



(13) 



((cr 2 )k(er 2 )k') = — — J d 3 pd 3 p' (cr p( 7 k _pO-p,OV_ p ,) 

= 25( 3 )(k + k') J d 3 pP a (p)P a (\k-p\) 



Taking V a to be scale- independent 16], 



27r Jl- 1 V" 



d 3 p 



kl 3 



(14) 



The subscript L~ l indicates that the integrand is set 
equal to zero in a sphere of radius L^ 1 around each sin- 
gularity, and the discussion makes sense only for L^ 1 <C 
k <C fcmax- In this regime one finds |l6| 



V a 2(k) = 4V 2 Jn(kL). 



(15) 



This expression gives the correct variance 
({a 2 — {a 2 }) ) = 2((T 2 ) 2 , confirming the consistency of 
the finite box approach. 

Now consider the bispectrum Bq = B a i . Repeating 
the calculation leading to Eq. (|14fl . we find eight terms 
which can brought into a common form by a redefinition 
of the integration variables, giving 



d 3 p 



p-ki| 3 |p + k 2 | 



(16) 



where again L _1 means that a sphere of radius is 
cut out around each singularity. The integral may be es- 
timated by adding the contributions of the three singu- 
larities, and comparing the result with Eq. ijSJ one finds 



20 V 3 

f NL = -A Y ^]n(kL). 



(17) 



with A = O(l) and k = min ki. Expressed in terms of the 
non-gaussian fraction r^ a = (V a 2 /Vq) 1 ^ 2 '', this becomes 
r Ca ~ 0.04/ NL 1/3 



3 



To compare with observations, we choose a minimal 
box size, so that ln(/cL) is roughly of order 1 on cosmo- 
logical scales. Since /nl has only mild scale-dependence, 
we can use observational bounds on that are obtained by 
taking /nl to be scale- independent. The current bound 
|/nl < 100 corresponds to r a i < 0.18. Planck [T^ will 
be able to detect an un-correlated Gaussian-squared com- 
ponent of the primordial curvature perturbation, pro- 
vided that its relative amplitude r a i is bigger than 0.04. 

The trispectrum. Computing the connected four-point 
correlator of a 2 we find 



d 3 P 



23 perms , 



IP - ki| 3 |p 



k 2 | 3 |p- 



4-24| 



(18) 



where the permutations are of {kik2k3k4} giving actu- 
ally 3 distinct terms. A sphere of radius L _1 is cut out 
around each singularity, and imposing k 3> L^ 1 on the 
wavenumbers in Eq. i|13[l gives |ky| ^ L -1 . The corre- 
spondence between Eqs. (I18fl and (|10fl is seen by drawing 
Feynman-like diagrams as for instance in 18j. Estimat- 
ing the integral by adding the contributions of the four 
singularities, and comparing with Eq. I|10|) . one finds 



tnl = 16S|§ ln(fcL) . 



(19) 



with B = O(l) and k = minjfci, |kj,„|}. For the fraction 
of non-gaussianity this corresponds to r a i ~ 10 _1 tnl 1 ^ 4 - 
This result might be competitive with the one coming 
from /nl, when tnl is bounded by WMAP. 

A Gaussian- squared isocurvature perturbation. In ad- 
dition to the primordial curvature perturbation, there 
might be a primordial isocurvature perturbation. Like 
the former, the latter is to be evaluated a few Hubble 
times before cosmological scales start to enter the hori- 
zon, and is specified by 



g Srii <5n 7 

Ti% Tizzy 



(20) 



where the n's are number densities and i runs over cold 
dark matter particles (C), baryonic matter (B) and neu- 
trinos (v). Observation requires roughly \Si\ < 0.1C- De- 
pending on its origin, Si(x) might be a constant multiple 
of C( x )j but it might instead be uncorrelated with it [22j. 
The classic example of the latter case is a CDM isocurva- 
ture perturbation caused by a perturbation in the axion 
field. In that situation Sc will be the sum of terms lin- 
ear and quadratic in the axion field [ly. Let us suppose 
that the latter dominates and that ( is Gaussian. On 
large angular scales, ignoring baryons, the temperature 
anisotropics are given in term of the primordial curvature 
and isocurvature perturbation through 



/ AT\ _ /AT\ f AT \ C 2 



\ T J 



sw 



V T 



Sc (21) 



The angle-averaged bispectrum defined e.g. in [9j will 
consist only of the isocurvature one i.e. Bi lt £ 2t £ 3 = 
Bg S< ~j^ t . Since the coefficients of £ and Sc are sim- 
ilar, the observational bound on Vsc/'Pq from non- 
Gaussianity will be similar to the one on V^/Vq^ and 
hence competitive with the existing bound. 

The logarithmic scale dependence. We have found that 
the spectrum, bispectrum and trispectrum all increase 
like ln(fcL) with k a representative wavenumber. This 
factor arises from the infrared divergence of Eq. JT3J, 
which in turn is due to the interference of very long wave- 
length components. We adopted a minimal box size mak- 
ing the factor of order one on cosmological scales, though 
it might still be significant on very small scales relevant 
for instance for primordial black hole formation. But ir- 
respective of practicalities, the appearance of the factor 
seems to contradict a basic tenet of physics concerning 
the use of Fourier series; that the box size should be ir- 
relevant if it is much bigger than the scale of interest. 

To see what is going on, we need to consider the generic 
perturbation which contains a linear term aa. The 
crucial point is that the mean (a) is supposed to vanish 
within the chosen box of size L. Suppose now that we go 
to box of size M <C L, still comfortably surrounding the 
region of observational interest. (It might be the whole 
observable Universe, or a much smaller region around 
us where say the distribution of primordial black holes 
has been detected.) Inside the small box the appropriate 
variable is a = a — (c) M , and we will denote expectation 
values inside the small box by {) M - Then 



G(x) - a (M)«r(x) + ((7 2 (x)-( ( 7 2 )) , 



(22) 



with a(M) = a + 2(a) M . Finally instead of locating 
the small box around a particular observed region, let 
its location vary so that (a) M becomes the original per- 
turbation a smoothed on the scale M . It is clear that 
the stochastic properties of the perturbations within a 
randomly-located small box will be the same, whether 
we use such a box directly or whether we calculate them 
using instead the big box. 

Let us see how this works for the spectrum, bispectrum 
and trispectrum of Q a . The spectrum is "P^ a = a 2 V a + 
V a 2 . Defined within a particular small box it becomes 

V^{M)^a 2 (M)V (y + V a 2. (23) 

Using Eq. (|15(l . the expectation value is 

(V C (M)} = AV 2 a (HL/M) + In(kL)) = V c . (24) 

These equations translate to the following words; going 
to the small box increases the the size of the cutoff around 
the singularities, from L _1 to M _1 , but this is compen- 
sated by the change in the coefficient a(M). 

Exactly the same thing happens for the bispectrum 
and trispectrum. Using our previous results, the bispec- 
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trum within the large box is 



B f 



,4„2t-,2 f 1 



a \ 1.31.3 
\ 1 2 



cyclic 



Jl- 1 



d 3 P 



L -ip 3 |ki-p| 3 |k 2 + p| 



(25) 



Replacing a by a(M) and L by M gives the bispectrum 
B(_(M) within a small box, leading to 



(B C (M)) = (2n) 3 Vl 



4tt ln(L/M) 
d 3 p 



1 



cyclic 



A/ 



-i p 3 |ki -p| 3 |k 2 + p| 



(26) 



The first term changes J M -i to J L -±, and a similar cal- 
culation works for the trispectrum. 

What we are finding here is a logarithmic 'running' of 
the (typical value of) the parameter a with the box size, 
similar to the running of field theory parameters with 
the renormalization scale. The logarithm would become 
a power if Va(k) went like a positive power. If more fields 
and/or higher terms in the expansion of £ were allowed, 
we would arrive at 'renormalization group' equations re- 
lating the running of the coefficients in the expansion. 
That is hardly of practical interest for the primordial 
curvature perturbation, but it will be applicable when 
the formalism developed in is extended to include 
non-Gaussianity. 

Returning to the primordial perturbation, consider 
again our assumption that £ (or Si) is quadratic in cr, 
corresponding to a ~ 0. If this assumption is valid within 
a box exponentially larger than the observable Universe, 
then it can be valid also for the minimal box that we 
adopted only if our location is to some extent untypi- 
cal. This situation was noticed for the axion [lil Ejf . 
providing an early example of a theoretically-motivated 
anthropic consideration. 

Conclusion. We have considered the effect of a sub- 
dominant primordial perturbation, which is uncorrelated 
with the main component. Unlike the main component, 
this one may be highly non-Gaussian, with the result that 
its stochastic properties are quite different from those of 
the main component. They have to be defined with re- 
spect to a comoving box of finite size, but a change in 
the box size is compensated by a change in the coefficient 
of the Gaussian component as measured by a typical ob- 
server, so that physical results are unchanged for such an 
observer. 

Coming to the observational side, and assuming that 
the sub-dominant is the square of a Gaussian, it will 
eventually be detectable if it contributes more than a 
few percent of the total. A smoking gun for this setup 
would be the detection of both primordial gravitational 
waves (indicating that the inflaton gives the dominant 
contribution to the primordial curvature perturbation) 



and non-Gaussianity (indicating that this contribution 
cannot be the only one). Another smoking gun would be 
provided by the observation of both the bispectrum and 
the trispectrum, the latter being at least equally impor- 
tant. 
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